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ABSTRACT 

We present a simple analysis of X-ray emission and absorption for hot-star winds, 
designed to explore the natural scalings of the observed X-ray luminosity with wind 
and stellar properties. We show that an exospheric approximation, in which all of 
the emission above the optical depth unity radius escapes the wind, reproduces very 
well the formal solution for radiation transport through a spherically symmetric wind. 
Using this approximation we find that the X-ray luminosity Lx scales naturally with 
the wind density parameter M/voo, obtaining Lx ~ (M/voo) 2 for optically thin winds, 
and Lx ~ (M /v 00 ) 1+s for optically thick winds with an X-ray filling factor that varies 
in radius as / ~ r s . These scalings with wind density contrast with the commonly 
inferred empirical scalings of X-ray luminosity Lx with bolometric luminosity LboI- 
The empirically derived linear scaling of Lx ~ ^Boi for thick winds can however be 
reproduced, through a delicate cancellation of emission and absorption, if one assumes 
modest radial fall-off in the X-ray filling factor (s ~ —0.25 or s ~ —0.4, depending 
on details of the secondary scaling of wind density with luminosity) . We also explore 
the nature of the X-ray spectral energy distribution in the context of this model, and 
find that the spectrum is divided into a soft, optically thick part and a hard, optically 
thin part. Finally, we conclude that the energy-dependent emissivity must have a 
high-energy cut-off, corresponding to the maximum shock energy, in order to reproduce 
the general trends seen in X-ray spectral energy distributions of hot stars. 

Subject headings: radiative transfer — shock waves — stars: early-type - 
stars:mass-loss — X-rays: stars 



1. Introduction 

One of the great surprises of the Einstein X-ray satellite mission was the observation that 
hot, luminous, OB stars are strong emitters of soft (~ 1 keV) X-rays (Harnden et al. 1979; Seward 
et al. 1979). Unlike the cooler, late-type stars, for which the X-ray emission was found to scale 
with the star's rotation, the observed X-rays from hot stars show a roughly linear proportionality 
to the stellar bolometric luminosity, L x ~ 10~ 7 Lb o i (Long &: White 1980; Pallavacini et al. 
1981; Chlebowski, Harnden, and Sciortino 1989). Observations from subsequent X-ray satellites, 
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most notably ROSAT, have generally confirmed this earlier scaling result, with some minor 
refinements. For example, from a ROSAT survey of 42 O stars, Kudritzki et al. (1996) find 
logLx/^Boi ~ —6.7 ± 0.35, with a somewhat tighter relationship when a weak dependence on a 
characterisic wind density M/vco is included, L x oc (M/foo) a38 L Bo i L34 - 

Even before these detections, X-ray emission from a narrow stellar "corona" was postulated 
to explain the superionization seen in the UV spectra of OB stars (Cassinelli & Olson 1979). This 
coronal model was further refined by Waldron (1984), but other studies have limited the extent, 
temperature, and fractional contribution of the total X-ray output of such purported O-star 
coronae (Cassinelli, Olson, & Stalio 1978; Nordsieck, Cassinelli, & Anderson 1981; Cassinelli & 
Swank 1983; Baade & Lucy 1987; MacFarlane et al. 1993). A much more generally favored 
scenario has been that these X-rays arise from shocks that form in the wind from strong intrinsic 
instabilities in the line-driving mechanism (Lucy & Solomon 1970; MacGregor, Hartmann, & 
Raymond 1979; Owocki k Rybicki 1984, 1985). Lucy & White (1980) and Lucy (1982) proposed 
phenomenological models in which wind material accelerated by this line-driven instability rams 
into ambient, "shadowed" material, forming a forward shock that accelerates and heats that 
material, producing X-rays. 

Dynamical simulation models of the nonlinear evolution of this instability have since indicated 
a quite different wind structure, dominated by reverse shocks that decelerate very rarefied, 
instability-accelerated material as it rams into compressed, dense shells ahead of it (Owocki, 
Castor, & Rybicki 1988). Since only a small fraction of wind material is actually heated by 
passage through such reverse shocks, the resulting direct X-ray emission is very low, generally well 
below the observed value (Feldmeier et al. 1997a). To overcome this limitation, subsequent wind 
structure simulations by Feldmeier et al. (1997b) introduced a chaotic, turbulent perturbation at 
the photospheric lower boundary. This seeds a wind structure with substantial velocity dispersion 
among the dense shells, which thus collide amongst themselves, yielding much more shock-heated 
material. For the specific case of a massive OB supergiant wind with rather strong and chaotic 
base perturbations, the X-ray emission produced in such models can be made to approach the 
observed values. However, the dynamical simulations for particular cases have thus far made no 
attempt to explain the broad scaling of observed X-ray emission with stellar luminosity. 

The goal of the present paper is to provide a firm basis for understanding these more general 
scaling properties of hot-star X-rays. Our approach here eschews detailed dynamical simulations, 
and instead examines phenomenalogical models of the X-ray emission, in order to define more 
clearly what overall properties are needed from instability-generated wind structure simulations 
to produce these observed scaling relations. A major point of our analysis is that a linear scaling 
of observed L x with does not arise naturally from expected properties of a wind-shock model, 
but requires a rather delicate balance of X-ray absorption and emission, which in turn requires 
a special form for the radial distribution of wind shocks. A further point is that such a balance 
between emission and absorption should be associated with specific trends in broad-band X-ray 
spectral properties. 
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The remainder of this paper is organized as follows. Beginning from the formal solution to 
the radiation transfer equation for X-rays emitted throughout an expanding stellar wind (§2.1), 
we develop a simple "exospheric" approximation through which we derive a general scaling law 
for X-ray luminosity from a simple power-law emission model (§2.2). By relating this to wind 
and stellar scaling laws we identify parameter requirements for reproducing the observationally 
inferred linear dependence of L x on L^i (§2.3). We next examine the energy dependence of the 
X-ray spectrum (§3.1) and then integrate the expression for the X-ray spectrum over energy in 
order to verify the result from the previous section that gave the linear Lx ~ Lb i scaling (§3.2) 
The final section (§4) summarizes our conclusions. In two appendices we examine the effect of 
assuming a constant wind velocity and verify the validity of the assumption of a smooth X-ray 
emissivity in terms of detailed numerical plasma emission models. 



2. Radiation Transfer Through an Expanding, Spherical Stellar Wind 

2.1. Formal Solution for Constant Velocity Expansion 

To develop a simplified description of X-ray emission and absorption in a spherically 
symmetric, expanding envelope, we initially consider monochromatic X-rays, deferring the 
discussion of energy dependence to the next section. Let us first define the usual (p, z) coordinates 
such that z is a coordinate along a ray, and p is the 'impact parameter', or minimum radial 
distance of the ray from the origin. From a formal solution of the transfer equation for X-rays 
emitted at a radius r = \Jp 2 + z 2 with emissivity i](r), the specific intensity observed at large radii 
(z — > oo) is 

I v = r r,{r)e~ TM dz. (1) 

J — oo 

The total X-ray luminosity is then given by integration over all rays 

roo 

L x = 8tt 2 I p pdp. (2) 
J o 

Here the ray optical depth is of the form 

roo 

T (p,z)= K x (r(p,z')) p(r(p,z'))dz', (3) 

J z 

where p is the mass density, and k x is the X-ray absorption cross section per unit mass. We note 
that the opacity arises in relatively cool, unshocked portions of the wind, and is due primarily to 
K-shell bound-free transitions in He, C, N, and O, Ne, Na, Mg, and Si. 

For an instability-generated shock model, X-ray emission is only expected once the wind has 
reached a substantial fraction (>50%) of the terminal speed Uqo, implying that the radial variation 
of mean wind density is within an order unity factor of that for a constant velocity (Owocki 1994; 
Feldmeier 1995). For simplicity, let us thus consider the case of expansion at a constant velocity 
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Voo (but see Appendix [A| for a treatment of a beta velocity law). Moreover, since ionization 
fractions are also fairly constant beyond about 0.5^00 (MacFarlane, Cohen, & Wang 1994), we 
assume also a radially constant X-ray mass-absorption coefficient, k x . The radial optical depth is 
then given simply by 

kJA f°°dr' R 1 
r (°> r ) = 7Z— / — = — > ( 4 ) 

where M is the mass loss rate, and 



47T-UOO Jr r' 2 



*.=^ (5) 

defines the radius of unit radial optical depth, i.e. r(0,i?i)=l. For nonradial (p ^ 0) rays, we 
obtain 



( \ p f°° dz> Rl 
rip, z) = R\ \ — 77 = — arccos 

Jz p 2 + z' z p 



—0(p,z), (6) 



\Jp 2 + z 2 

where is the angle between the +z ray direction and an outward radius vector. This suggests it 
would be convenient to recast the intensity integral (|]) in terms of this angle, 

V(r(p, 6)) 
/o sin 2 6 e 



PL „: n 2g p fl l0/p d0 - ( ? ) 



To proceed, we next note that X-ray emission arises primarily from two-body processes - 
recombination and collisional excitation - and thus has a volume emissivity rj (ergs s _1 cm -3 ster -1 ) 
that is proportional to the square of the density, 

ri(r) = iC s f m p 2 {r) (8a) 
= iC 2 J vP 2 {r), (8b) 

where f m is the mass fraction, or mass filling factor, of ambient wind that is X-ray emitting, and 
C s =p s /p is a factor to correct for any differences in the density p s of shocked material relative to 
the ambient wind density p (which assumes a smooth wind, given by p=M/47rr 2 f). We note that 
the factor C s accounts both for the post-shock compression and the deviations of the pre-shock 
density from the ambient wind density. The mass filling factor, f m , is related to the volume filling 
factor by f v = f m /C s . The emission coefficient £ (ergs cm 3 s _1 g -2 ) is related to the commonly 
used cooling function A (ergs s _1 cm 3 ) (e.g., Raymond & Smith 1977) by £ = A/Airp e p p , where 
p e and p p are the mean mass per electron and proton, respectively. Under typical stellar wind 
conditions £ is independent of density (at least when considered at the coarse spectral resolution 
associated with current X-ray detectors), and depends mainly on the electron temperature in the 
post-shock region. Assuming for now that the factor £C 2 /^ is spatially constant, the intensity 
integral (|7j) can be rewritten and evaluated as 



2(C - f '% ( 1 r{+^ 1 ' (9b) 
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Application to the luminosity integral (0) then yields 



a f°° 1 - e- 2lTX 



* g.isn'ic'iu 1 ^ (10b) 
= 2.44.^^, (10c) 

^X V(yQ 

where the approximate equality follows from direct numerical evaluation of the integral. Note 
that, although the wind volume emission varies as the density squared, the X-ray luminosity 
escaping the wind scales only linearly with the density parameter M/voq, due to the effect of wind 
attenuation. 



2.2. Exospheric Approximation for Power-Law Radial Emission 

It is of interest to compare this numerical result with the heuristic formula from an "exospheric 
approximation" (Cohen et al. 1996) that estimates the X-ray luminosity from the volume integral 
of the outward (i.e. radiated into 2ir ster) emission beyond the radius with unit radial optical 
depth, 

L x « Svr 2 r v (r) r 2 dr = 2?r ^ M (n) 

J Ri KxVoo 

where the latter equality follows from straightforward analytic integration. The rough agreement 
with the numerical result (eq. fl~0[] ) suggests that this simple exospheric formula should be quite 
useful for estimating L x in more complicated cases for which the required integrals are difficult to 
evaluate. 

Let us thus next estimate L x for the somewhat more general case when the factor £,C 2 f v 
varies as a power-law beyond some lower boundary radius R Q for X-ray emission, 

n(r) = p 2 (CC 2 J V ) (-^y ;r>Ro (12a) 

= ; r<R D , (12b) 

where (£,C 2 f v ) is the emission factor at r = R Q . Unlike for the special case R a = R* and s = 
considered above, exact evaluation for L x would now generally require numerical evaluation of a 
double integral over p and z. However, application of the exospheric formula (eq. |IT|) readily 
yields an approximate analytic scaling law for the X-ray luminosity, 

L x « (—) 2 R -%m a x[R ,R 1 ]y-\ (13) 

2(1 -s) yvoo J 

We note parenthetically that this can be recast as a fitting formula, 

Lx - ( 4^EM W ( % ) 13 , (14) 
2(1 — s) \m&x(R , R\) J 
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where EM W = J^n 2 dVol is the total wind emission measure above R Q . Alternatively, we can 
express this in terms of the basic stellar wind parameters, 



^ « ( —) S I — 1 ; ^ < i?i (optically thick) (15a) 



; R > Ri (optically thin). (15b) 



2(1 - s).R V it. 
(&f v ) (m\ 2 



2(1 -s)R \v c 



Equations ( |15b| ) and ( |15b ) represent the X-ray luminosity for, respectively, the 'thick' and 'thin' 
wind cases. Note that the luminosity now scales with (M/i^) 2 for the thin wind, but only with 
(M/uoo) 1+s for the thick wind. 

We reiterate that these results have all been derived for monochromatic X-rays. We will see in 
§3 that a given stellar wind can be optically thick and some energies and optically thin at others, 
leading to spectral structure based on the energy-dependence of the X-ray attenuation. 



2.3. Relation to Wind Scaling Laws 

The above simple analysis illustrates that the most direct scaling of the X-ray luminosity 
should be with wind density, and not, as has generally been inferred empirically, with the 
bolometric luminosity. Let us thus next consider whether this inferred scaling could be explained 
through secondary scaling of wind density with luminosity. Within the standard Castor, Abbott, 
and Klein (1975; hereafter CAK) theory for line-driven winds, the mass loss is predicted to scale as 

M ~ L Bo i 1/Q ' Mlj}' a ', (16) 

where M e ff is the effective stellar mass {M e ff = M*(l — TEdd))i an d oc' = a — 8 is a combination of 
the two line parameters, a and 5 defined by CAK and Abbott (1982), which respectively describe 
the distribution of line opacity and its dependence on ionization level. 

As first noted by Kudritzki, Lennon, and Puis (1995), this theoretical scaling is in good 
general accord with an empirically derived 'wind-momentum-lumosity' relation 

MvooRl 12 ~ L Bo i 1/Q ', (17) 



wherein the connection with the CAK mass loss formula (16) follows from the scaling of the wind 
terminal speed with the surface escape speed, 



Meff/R*, (18) 



which is a well known theoretical (CAK; Pauldrach, Puis, & Kudritzki 1986) and observational 
(e.g., Lamers, Snow, & Lindholm 1995) result. For O supergiants, the extensive analysis by Puis 
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et al. (1996) empirically obtained a wind- momentum-luminosity relation with a slope implying 
cJ = 0.57, which thus gives 

M L B oi L75 i# 5 



M*{l-T E dd) 



(19) 



The dependence on the Eddington factor is significant only for stars with very high 
luminosities. Between an early-B star like A Sco and an O supergiant like £ Pup, (1 — VEdd) 
changes by only a factor of 2, while the bolometric luminosities vary by almost three orders of 
magnitude. 

To convert eq. (^) to a scaling with luminosity alone, we can ignore the weak dependence on 
stellar radius; but we may want to allow for a systematic trend between luminosity and mass. For 
example, from stellar structure theory, one expects a quite strong trend e.g. L-q q \ ~ M^- 5 , yielding 
M/v oo ~ -^Bol ^ ■ The observed Lx ~ -^Boi relation would then be reproduced by the optically 
thick wind scaling relation ( |15b| ) if s ~ —0.25. 

By contrast, if we simply ignore any systematic trend between mass and luminosity, the 
observed Lx ~ ^Boi relation would be reproduced by a filling-factor index s « —0.40. Either 
scaling is qualitatively consistent with detailed numerical simulations of the line-driven instability, 
which generally show a gradual drop-off in X-ray production with radius beyond the point where 
strong shocks begin to form in the wind. But overall it seems that, with such a wind-based model 
for X-ray emission, connection with the observed Lx ~ -^Boi is inherently indirect, requring a 
rather specialized cancellation between the wind emission and absorption. 

Finally, we note that for lower-density winds that are optically thin to X-rays, the X-ray 
luminosity becomes independent of s, scaling as 

2-7 _ t 3.5 




~W ' or L B oi , (20) 

where the two scalings depend on whether one includes or ignores the systematic variation of 
luminosity with mass. Either form is in general accord with empirically inferred X-ray scalings, 
which show the dependence of X-rays on bolometric luminosity becomes much steeper around 
early-B stars, corresponding roughly to where the winds are becoming optically thin (Cohen et al. 
1997). 



3. Energy Dependence 

3.1. Power-law Absorption and Emission 

All of the above analysis is for monochromatic X-rays. However, current instruments 
{Einstein, ROSAT, ASCA) are sensitive to X-rays over almost two orders of magnitude in energy, 
and even with their quite modest spectral resolutions, have provided some basic information 
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about the spectral energy distributions of X-rays in hot stars. So let us now consider the energy 
dependence of the absorption and emission, and of the resulting luminosity spectrum. Given the 
modest energy resolution of X-ray observations to date, we focus here on the broad-band spectral 
properties. The X-ray absorption - apart from several prominent K-shell edges - can be fit roughly 
by a power-law in energy, 

^(E x ) - ||. (21) 

Hillier et al. (1993) find, for £ Pup (04 If), a ~ 2 in layers where helium remains doubly ionized, 
and a ~ 2.8 where it has recombined to He + . For e CMa (B2 II), detailed calculations give a ~ 1.8 
between 100 eV and 1 keV (Cohen et al. 1996). This energy dependence of opacity means a given 
wind can be optically thick at some energies, and optically thin at others. By setting Ri = R Q , we 
can derive a critical energy, 

that separates the parts of the spectrum for which the wind is optically thick (E x < E\) from 
those for which the wind is optically thin {E x > E\). 

To derive the form of the luminosity spectrum one must specify the energy dependence of 
emission as well as the absorption. For simplicity let us again assume that the emissivity has a 
power-law energy dependence, 

^E x ) = |-, (23) 



which means that equation (12a) can be written as 77 = p 2 (^Cgf v ) E~ b (r/R ) s . In Appendix || 
we show a comparison of this simple, smooth emissivity model to detailed numerical simulations 
of collisional equilibrium plasmas. This comparison shows that at the low resolution of current 
instruments, the line-dominated spectra are surprisingly well approximated by a smooth function 
such as equation (|23|), with a power-law spectrum being produced by a power-law temperature 



distribution. The luminosity spectrum is then given by, 



at, f E \ a ( 1 ~ s )- b 

M=H ;E X <E X (24a) 



dE x V E\ 

= Lx b ; E X >E U (24b) 

where the luminosity at the energy at which wind optical depth is unity, L\=L X (E\), has the 
form 

L = (^f v ) (*kB*\V* (M\ 2 - b/a 
1 2R {l-s) \ k xo ) \ Voo ) ' 1 ' 

While equations ( |24a| ) and ( |24b| ) are convenient for plotting the spectrum for a given model, 
we can alternatively write this in a form that more explicitly displays the scalings with stellar 
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parameters, 

i-» / w\ 1+ 



2 




(C^ s 2 /,)o fM , ^ ; E x >E l7 (26b) 



2i? (l-s) U c 



where, again, equations ( |26a| ) and ( |26b| ) represent the expressions for, respectively, the optically 
thick and optically thin wind domains (cf. eqs. [|15b|| and ||15b||). 



Using the simple power-law absorption and emission models, the energy-dependencies in the 
two domains have slopes that differ by a(l — s). Furthermore, if a(l — s) — b > then E x = E\ 
is also the energy of peak flux, and the value of this peak is dictated solely by the absorption 



properties of the wind. Given equation (22) one would expect the observed values of E\ to vary 
significantly from low-density B star winds to high-density O star winds. However, although 
some variation in the peak intensities are seen, it is much less dramatic than this model predicts. 
In addition, the relative slopes of the optically thin and optically thick parts of observed X-ray 
spectra are generally seen to differ by much more than a(l — s), given reasonable values of a and 
s. These trends are shown in Figure |l] for a representative O star (£ Pup) and a representative 
early-type B star (A Sco). 

It is thus apparent that a universal power-law emissivity function is too simple an assumption. 
An obvious problem with this functional description is that it implies the presence of plasma 
having an arbitrarily high temperature. It therefore seems appropriate to consider a truncation to 
the power-law emissivity introduced in equation (£3) that corresponds roughly to the maximum 
shock energy in the wind. We therefore modify the emissivity as 

£{E X ) = ^ e - E */ E ™*, (27) 

where E max is the maximum plasma temperature, corresponding to the maximum shock energy. 
The luminosity spectrum is now given by 

AT / J? \ 0,(1— s)— b 

dL x L (^\ e -E x /E max . Ex < Ei (2ga) 



dE x \EiJ 

= Ll (j^j e -E*/Ema* . Ex >E U (28b) 

where L\ has the same definition as before. The exponential cut-off in this model brings the 
luminosity spectrum more in line with the data shown in Figure |l], as well as that for most other 
OB stars. 

With this general expression for the X-ray spectra, one sees that there are two different types 
of spectra that can occur. If E max > E\ then the high-energy fall off of the spectrum occurs 
above E max , where the spectrum would be decreasing as a power-law already. Relative to the 
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Fig. 1. — BBXRT (detector A) spectrum of £ Pup (04 If) taken from Corcoran et al. (1993) (a). 
This spectrum has been "unfolded" so as to correct for the energy-dependent sensitivity of the 
detector. Note the peak flux is near 0.7 keV, once the oxygen K-edge feature near this energy is 
corrected for. Given the parameters of £ Pup (M = 5 x 10 6 M yr _1 ; = 2.2 x 10 3 km s _1 , 
R* = 1.3 x 10 12 cm) we find a value of E\ ~ 2 keV, however. Note also that the power-law indices 
of the slopes of the low-energy and high-energy portions of the spectrum differ by at least six. The 
ROSAT PSPC spectrum of A Sco (B1.5 IV) (b) has a peak energy (corresponding to E±) that is in 
much better agreement with the predicted value of 0.1 keV. The slightly higher observed peak value 
is entirely consistent with the finding that the mass-loss rates of B stars have been systematically 
underestimated (Cohen et al. 1997). 
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simple power-law expression model, this spectrum has a peak determined by E%. Effectively, then, 
the main change is to cause a much steeper high-energy fall-off in the spectrum. However, if 
Emax < El then the peak in the spectrum is at a softer energy than it otherwise would have been 
(at an energy dictated by E max rather than E\), and the fall-off in the spectrum with increasing 
photon energy occurs that much sooner. 

When O star X-ray data are fit with two-temperature models (e.g. Hillier et al. 1993), the 
higher temperature component tends to be approximately 5 x 10 6 K, rarely exceeding 10 7 K. Even 
for the cooler B stars, the highest temperature component derived from model fitting is roughly 
5 x 10 6 K (Cassinelli et al. 1994), indicating that E max should be relatively constant for hot stars, 
with a value near 0.5 keV. This is in contrast to the values of E±, which vary by more than an 
order of magnitude from mid-B to early-0 stars. Therefore, which of these two cases applies to a 
given star depends primarily on the value of E\, which is governed largely by the mass- loss rate. 
The case with E max < E\ is thus appropriate to O stars, with their optically thick winds, while 
the case with E max > E\ is relevant for B stars, with their thinner winds. For these low-density 
winds, the peak energy in the spectrum is determined by the value of E\, whereas for the higher 
density, optically thick winds the peak of the spectrum is relatively constant at the universal value 
of E max « 0.5 keV. 



3.2. The Lx/LboI Scaling Relation from the Full Spectrum 

With the parameterization of the energy-dependence of wind attenuation and wind emission 
introduced in the previous subsection, we can return to the explanation for the scaling relationship 



discussed in §2^3. Because the O stars and early-B stars for which the Lx/L-q q \ scaling holds have 
optically thick winds with E max < E\, the optically thin portion of the wind (E x > E±) contributes 
negligibly to the overall X-ray luminosity and we therefore can consider only the optically thick 
expression, given by equation fl28a|) . Furthermore, because of the negligible contribution of the 
high-energy portion of the spectrum, we can carry out the integration over all energies without 
significantly affecting the result. The total X-ray luminosity emergent from a hot star wind in this 
model is then given by 

This expression preserves the Lx ~ (M/i>oo) 1+s dependence which yields Lx ~ -^Boi with 
s ~ —0.25 or s ~ —0.4, as we derived in § |2.3] for monochromatic X-rays. We note that this analysis 
holds only for the optically thick winds that correspond to O stars which have Lx/^Boi ~ 10~ 7 . 

The relationship for monochromatic optically thin winds, Lx ~ (M/foo) 2 , is also consistent 
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with the energy-dependent treatment, so long as E\ is below the soft edge of the instrumental 
bandpass. Otherwise the observed spectrum has significant contributions from both an optically 
thin power-law region and an optically thick power-law region. In this hybrid case, the dependence 
of the X-ray luminosity on the bolometric luminosity is complicated and falls somewhere between 
the linear dependence of the thick winds and the much steeper dependence of the thin winds. 
This would cause a less discontinuous fall-off in the Lx/^Boi relationship than is inferred from the 
monochromatic analysis in §2.2. 



4. Conclusions 

Despite its approximate nature, the analysis described here provides a useful initial framework 
for examining the interplay among the basic processes of X-ray emission and absorption that 
determine the observed X-ray spectra from a hot-star wind. 

Some of the main points are: 

1. The simple procedure (cf. eq. fl"lf ) of including only the X-ray emission above the radius 
R\ of unit optical depth (the exospheric approximation) provides an accurate means to 
account for X-ray attenuation effects. A key aspect of this approximation is the assumption 
that v ~ Voo, which is indicated theoretically. If the X-rays are formed deeper in the wind, 
then this constant-velocity approximation will underestimate the X-ray luminosity, but in 
Appendix |A| we estimate that the error is generally no more than about 50%. 

2. For optically thin winds, the X-ray luminosity increases as (M/voo) 2 , as would be expected 
based on the density-squared sensitivity of the thermal emission processes. However, the 
dependence of the X-ray luminosity on wind density is less steep in the case of optically 
thick winds, with L x ~ (M /i>oo) 1+s , where the X-ray filling factor has a power-law radial 
dependence, / ~ r s . 

3. Indeed, the natural scaling of the x-ray luminosity in hot-star winds is with the wind density 
parameter M/voo, and not with stellar parameters such as the bolometric luminosity. 

4. However, if one introduces a radial power-law scaling of the filling factor, then the observed 
Lx ~ LboI relation can be recovered with a value of s ~ —0.25 to s ~ —0.4 for the radial 
power-law index. 

5. The energy dependence of the unit-optical-depth radius R± can play an important role in the 
shape of the energy spectrum. Indeed, in the simple shock model presented here, the peak of 
the spectrum occurs at the energy E\ for which X-rays become optically thin at the initial 
shock onset radius R Q , as long as the emissivity does not begin to fall-off rapidly at energies 
below E\. This condition should be satisfied for the low-density B star winds. However, for 
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O star winds, the peak of the spectrum is governed by the high-energy cut-off in the shock 
strength distribution. 

6. The X-rays observed by satellite telescopes are, especially for stars with optically thick winds, 
only a fraction of the total X-ray production. A significant amount of the X-ray production 
is absorbed by the optically thick winds of these stars. In addition, a significant amount 
of the total generated X-ray emission may fall outside of the observational bandpasses, 
especially in the very soft X-ray and EUV. 

Future work will apply the principles illuminated here to interpretation of more fundamental, 
numerical simulations of X-ray emission from instability-generated shocks, to a detailed analysis 
of shock heating and cooling mechanisms, and to the interpretation of observational data. This 
will help clarify the strengths and weaknesses of the wind instability paradigm for explaining the 
observed X-ray spectra and scaling relations. 

This research was supported in part by NASA grant NAG5-3530 to the Bartol Research 
Institute at the University of Delaware. 



Although there is good reason to believe that the X-ray production in hot star winds occurs 
primarily in regions of the wind where the velocity is close to the terminal velocity, it is possible 
to relax this assumption and still analytically solve for the X-ray luminosity under the exospheric 



where the emissivity r\ ~ p 2 , and R\ is the radius at which the radial optical depth is unity, 
r(0, = 1. For a 'beta' velocity law of the form 



A. Exospheric Model with a Beta Velocity Law 



approximation (eq. [0]), if the filling factor is radially constant. 



The exospheric approximation is given by 




(Al) 




(A2) 



the radial optical depth is given by the integral 




(A3) 



which through a variable substitution y = (1 — R*/r) can be readily evaluated as 




(A4) 
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Setting r = 1 leads to a new expression for the unit-optical-depth radius, 
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where R1=k x M/4ttv oo is the previous result for a constant velocity flow, which simply represents 
the special case (3 = 0. 

Making the same change of variables in the exospheric expression for the X-ray luminosity 
and using the newly defined lower limit of integration, R\ eta , one finds 
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Scaled by the constant- velocity exospheric expression (eq. |ll|), this can be cast as a correction 
factor, 
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Note that as R®/R* approaches infinity, L x eta /Lx approaches unity. Furthermore, for (3 = 0.8 
if Ri = 2R* then L b x eta /L x » 1.25. Choosing an E*[ that gives R\ eta = 1.5i?*, which seems a 
reasonable minimum radius for onset of instability-generated X-rays, L b x eta is just over 50% greater 
than Lx- 



B. Comparison of Power- law Emission Spectra to Numerical Models 

Given that the optically thin, thermal emission expected from the shock-heated regions of 
hot-star winds is dominated by lines (e.g. Raymond & Smith 1977; Mewe, Groenenschild, & van 
den Oord 1985), we may ask how well the smooth power-law models we have assumed in this paper 
represent the true spectra. We thus created a series of collisional equilibrium spectral models 
using the MeKaL code (Mewe, Kaastra, Liedahl 1995) and assuming B-star abundances (Geis & 
Lambert 1992). These are power-law differential emission models of the form ~ 1/T C that 
are constructed from a superposition of multiple isothermal models. The individual isothermal 
models used to make the power-law differential emission models have temperatures on the interval 
5.0 < logT(K) < 7.7. The spectra were calculated on the range 15 eV to 2500 eV. Although 100 
eV is generally taken to be the lower bound of the X-ray bandpass, significant emission occurs at 
lower energies, even for very hot plasma. We have compared these single power-law temperature 
distribution models to the single power-law spectral models described by equation (| 



Currently, the best spectral resolution available from space-based observatories is A/ A A ~ 20 
from BBXRT and ASCA. We therefore binned the detailed models at slightly better than this 



-15- 



resolution, and fit the resulting spectra with single power-laws, £ ~ E~ . Three examples of this 
exercise are shown in Figure ^. The power-law assumption is seen to be quite good. Reasonable 
X 2 values result from these fits if statistical uncertainties of only a factor of 2 to 3 are assumed for 
each bin. It should be kept in mind that the actual resolution of most available X-ray data sets 
is significantly worse than the binning of these models indicate, so that the observed data are, in 
fact, even more well-described by power-laws. 

There are several interesting things to note from this exercise. First, the global structure in 
these model spectra are dominated by the emission bump centered near 0.8 keV. This feature 
is due to thousands of iron L-shell emission lines. Second, a very significant fraction of the flux 
emerges not in the soft X-ray, but actually in the EUV, below 100 eV. Although this radiation 
is very significant (not least because of the role it plays in determining the wind excitation and 
ionization conditions), it has been directly observed in only one hot star, e CMa (Cohen et al. 
1996). Finally, the best-fit spectral power-law slopes are generally steeper than the power-law 
slopes of the differential emission measure (plasma temperature) distributions by roughly 0.7. This 
factor is a reflection of the fact that on the temperature range 10 5 < T < 10 7 , hotter plasma emits 
less efficiently than cooler plasma. Indeed, the slope of the frequency-integrated cooling-curve 
between the temperatures of 10 5 and 10 7 K (e.g. Raymond, Cox, & Smith 1976) is very close to 
the value 0.7. 



10 a 

hi/ (eV) 



10 3 



Fig. 2. — Temperature power-law spectral models (solid lines) calculated from the MeKaL code are 
shown along with the power-law spectral fits (dotted lines) . The best-fit power-law spectral indices 
(b) are indicated along with the differential emission measure slope used in making the models (c). 
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